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Abstract 

We demonstrate the equivalence of a Non-Markovian evolution equation with a linear memory- 
coupling and a Fokker-Planck equation (FPE). In case the feedback term offers a direct and 
permanent coupling of the current probability density to an initial distribution, the corresponding 
FPE offers a non-trivial drift term depending itself on the diffusion parameter. As the consequence 
the deterministic part of the underlying Langevin equation is likewise determined by the noise 
strength of the stochastic part. This memory induced stochastic behavior is discussed for different 
initial distributions. The analytical calculations are supported by numerical results. 
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I. INTRODUCTION 

One century after Einstein's explanation of Brownian motion compare also j^, 
stochastic processes are ubiquitous in almost every branch of physics |^ and beyond it 
like in economics ^ , chemistry ^ and biology |6] . The notion " Brownian motion" not only 
stands for the movement of pollens suspended in water, but is a generic term for the Wiener 
p..oces. ,7, see also the .ex.booUs gfa. Two .a,, approaches de.cnb,.g s.ochas* 
processes has been developed, related to the Langevin and to the Fokker-Planck equation, 
respectively. The underlying basic for the stochastic description is the ability of a separation 
of time scales leading to an evolution equation with a deterministic part supplemented by a 
stochastic force. In this Langevin equation the probability of the stochastic force is in general 
an input. An alternative way is the treatment of the stochastic process by the Fokker-Planck 
equation (FPE) 0, 0, Q , where a parabolic partial differential equation for a probability 
density is the basis of calculation. As already known both description are equivalent for 
special stochastic forces 0, [l^ [ill . In the last decades there is an increasing effort in gen- 
eralization of the L angevi n equation respectively the Fokker-Planck equation by inclu ding 



memory effects I^. ll3lll4l ll5ll . In the framework of investigation of anomalous diffusion [H 
and fractional diffusion [13] the descriptive equations are non-Markovian (non-local in time). 
The combination of diffusion with feedback couplings has not only an influence on the long 
time behavior, for the stationary state, but on the dynamics in an intermediate time regime. 
Such a memory dominated behavior is well established in analyzing the freezing processes 



in undercooled liquids 
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iali8, 



19|, where the underlying mathematical representation is 
based on a projector formalism proposed by Mori 0. In Q two of us studied a simple 
model, which still includes the dynamical features of evolution models as conservation of the 
relevant quantity p{f, t), which could be interpreted as probability density for a particle and 
moreover, a time delayed feedback coupling. Due to this coupling the generic behavior may 
be changed by the additional delay effects and could lead to non-stationary steady state 
solutions. 

The aim of the present paper is to show the equivalence between a special non-Markovian 
equation j2l| and the standard FPE with a non-trivial drift term. Firstly we present the 
non-Markovian model and its main features which consists of the occurrence of a non- 
trivial stationary solution and the dependence of the solution on the initial distribution 
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Po{r) = p{f^ t = 0). As the result the underlying FPE reveals a drift term which is itself de- 
termined by the strength of the stochastic force. The corresponding potential are discussed 
in detail. The analytical results are supported by numerical simulations. 



II. THE NON-MARKOVIAN FOKKER-PLANCK EQUATION 

In media with a spatial-temporal accumulation process transport phenomena should be 



described by a stochastic approach [2^ based on probabilities. The time evolution of the 
probability density could depend on the history of the sample to which it belongs, i. e. the 
changing rate of the probability should be influenced by the changing rate in the past and 
so the evolution equation of the probability has to be supplemented by memory terms. A 
recent overview is given in ^Sj- The obvious modification of the FPE by memory effects is 



to replace the conventional equation by j21| 

t oo 

dtp{f,t) = M{f,t;p,Wp) + J dt' j d'^r'IC{f-r,t-t';p,Vp)C{r,t';p,\/p). (1) 

~oo 

This equation is of convolution type and consists of two competing parts standing for pro- 
cesses on different timescales. The first part manifested by the operator Ai characterizes 
the instantaneous and local process, whereas the second part with the operators /C and C 
represent the delayed processes, the memory. In general all the operators may be non-linear 
in p{f, t) and Vp(r, t). Physically it means that the time scale of the memory is determined 
by the relevant probability p itself. The specification of the operators has to be according 
to the physical situation, which one deals with. However the main feature of the probability 
density p(r, t) is its conservation in time: 

dP{t) d 



oo 



dt dt 



d\p{f,t)=0. (2) 



To preserve p the instantaneous term Ad has to be related to a probability current, e. g., 
Ai oc V-j. For an arbitrary polynomial kernel K{f, t) the conservation law (0) is not fulfilled 
in general. Using Laplace transformation one can show directly, compare for details 2l| , that 
the choice C = — dtp{f, t) guarantees conservation. Thus our starting equation is written in 
the form 

t oo 

dtp{f, t) = D V^p{f, t)- j dt' j d'^r'K{r- r*, t - t') dfpi'F', t') (3) 



as a special realization of Eq. (^. This Fokker-Planck equation relates p{f,t) to p{f^t') 
with < t' < t unlike to a conventional one, where the evolution is only dependent on the 
probability at present time. Moreover (jH)) offers a coupling between dtp{f,t) and dt'p{f,t'). 
The mixing of time scales leads to a substantial modification of the long time limit. To 
demonstrate the modification let us consider a very simple choice by a strictly spatial local, 
but time independent kernel 

K{r,t)=fiS{f), (4) 

where parameter /i > characterizes the strength of the memory. By this choice the spatial 
and temporal variables are decoupled. Inserting the kernel in Q one gets 

dtp{f, t) = D V'^p(f, t) — fx [p{r, t) — Po{r)] with Po{r) = p(r, t = 0) . (5) 

The time independence of the kernel means that all times t' (0 < t' < t) in the past have the 
same weight and so there is a very strong memory with a direct coupling of the instantaneous 
value to the initial value. The memory induced feedback to the initial value appears as a 
driving force. Without this coupling one can interpret the equation as a description of 
a particle, which performs a diffusive motion, where the probability density p(r, t) decays 
on a time scale . As (0) is a linear equation and so the solution of it could be found 
analytically for arbitrary initial conditions 



(6) 



-oo L 

where G{r,t) is the Green's function of the conventional diffusion equation 

From the general solution, some properties could be followed easily such as if the initial 
distribution is non-negative poif), so the p{r,t) does provided /x > 0. The second moment 
s(t) could be calculated 



oo 



s{t)= j r'p{r,t) d'^r = '^^^^^ j d'^r + j p^{r) df^r . (7) 

— oo — oo — oo 

Notice that for the limit of vanishing memory one can easily verify that the last 

equation shows conventional diffusive behavior. The selection of the initial distribution is 



the essential point in our model and so three example are given to illustrate the solution of 
(0). Without lack of generality we concentrate our calculation on the one-dimensional case. 
It can be generalized to higher dimensions. 



III. RESULTS 

Obviously the results for a system with memory should be sensitive with respect to the 
initial distribution or at least from configurations in the past. Therefore, we study different 
reahzations for po{r) separately. 

A. Stationary solution 

The first example is the delta-starting distribution po{x) = poS{x). Substituting this in 
(P) the following solution is calculated 

t) = -jp= A^'"-^) + ^ [U{x- D, /x) + f4x; /.)] (8) 



itK X 



erf I , + ^/JIt \ — sgn(±x) 
\VlDt J 



(9) 



where erf(x) is the error function j23| and k, = The first part is the solution of the 

homogeneous equation showing temporal decay with time constant yU~^. In the long time 
limit the system shows a non-trivial stationary solution 

lim t) = ps{x) = ^ e-"l"l . (10) 

Such a stationary solution is due to the permanent coupling to the initial distribution and 
the greater yU > the stronger is this effect and more pronounced are the deviations from 
the pure diffusive behavior (/i = 0). The generalization to higher dimensions < < 5 and 
arbitrary initial condition can be directly calculated. It results in 

p,{r) = ^ / rfV P'^' (k If - f I) , (11) 

(27r)2^ \f-r'\~ ' 

with = fi/D and Kt,{x) is a modified Bessel function 
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Using the last result in Eq. (fTT|l we get for the Gaussian distribution po{x) = Po^ 
the following result 



Ps[X) 



4 V A 



(12) 



with 



g±{x; P, A, k) = e=^'^^ erfc ± xy/X^ and /5 



(13) 



4AZ) 2v^ 

Here erfc(x) is the complementary error function 23]. In case of an exponential decreasing 
initial distribution po{x) = poe~^'^' the calculation leads to 

~ A I X I 1"^ I " 



Psi\x\) = < 



Po ■;;t^ 



for A 7^ K 



for \ = K 



(14) 



Let us note that the stationary distribution depends on all cases on the initial conditions. 



B. Relationship to Fokker-Planck Equation 

The conventional form of the FPE including an external force, has in the one-dimensional 
case the following form 

where D is the diffusion coefficient, supposed to be constant here, and /(x) is the drift force, 
for which f{x) = —dU{x)/dx with U{x) as corresponding potential. On the one hand, the 
diffusion coefficient D measures the intensity of the noise and represents the stochastic part of 
motion, whereas the drift coefficient f{x) corresponds to the force experienced by the system 
and so it describes the deterministic part of motion. In this subsection we calculate the force 
f{x) which corresponds to the potential U{x) for the different starting distributions Po{x), 
in such a way that both Fokker-Planck equations (0) and ()15|) are equivalent. To do this 
deterministic parts of both equations are compared in the long time limit, in the stationary 
state. Provided Ps{x) ^ the formal solution is found by integration 

/i / [Psii) - Poii)] di 

/(^) = ^ n + ^ , (16) 
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(a) Six) 



(b) exp {-x'^) 



(c) exp(- \x\) 



FIG. 1: Starting distribution po{x), po = 1 and A = 1. 



where C is an integration constant, which one could set to zero. To show the equivalence 
mathematically correct, one has to do this comparison in the following way. First take 
an arbitrary function h{x) with bounded support, then integrate the product of h and the 
deterministic part of (0) resp. (fT3j) over the complete real line, and finally compare the 
results of these integrations. If both integrations are equal, then the functions are equal. In 
case of arbitrary dimension the a- component of the drift force is given by 



Here the integration constant is assumed to be zero as in the one dimensional case. Fi- 
nally we present the results for three different initial distributions, depicted in Fig. ^ the 
corresponding stationary solutions is shown in Fig. |2l the drift term and the corresponding 
potential in Fig. Eland Fig. HJ respectively. For the delta-like starting distribution the drift 
term can be calculated to 





(17) 



and so one can verify the following potential 



u{x) = ^/JId I X 



(18) 



In case of a Gaussian initial distribution an analogue calculation leads to 




(19) 



whereas an exponential starting distributions yields to the drift term 




for K X 



for K = X 



(20) 
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(a) ^(a;) (b) exp {-x^) (c) exp(- \x\) 

FIG. 2: Stationary state Ps{x), po = 1, X = 1, D = 1 and /x = 2 (units arbitrary). 




(a) 6{x) (b) exp (— a;^) (c) exp(— |a;|) 

FIG. 3: Drift coefficient f{x), po = 1, X = 1, D = 1 and n = 2 (units arbitrary). 

The underlying potentials can be obtain after an integration. The results are shown in Fig. 4. 
It should be noticed that the potential is likewise determined by the diffusion constant D 
which is a measure of the stochastic force. This point will be discussed in the next section. 




(a) ^(a;) (b) exp (— a;^) (c) exp(— |a;|) 

FIG. 4: Potential U{x), po = 1, X = 1, D = 1 and n = 2 (units arbitrary). 
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C. Deterministic motion in the potential 



As discussed before the speciality of the present approach consists of including the dif- 
fusive parameter D likewise into the force term f{x) of the FPE. Insofar it seems to be 
reasonable to study the corresponding Langevin equation, especially in the small noise limit. 
The trajectories obey the equation 
dx{t) 



dt 



f[x{t)]+i{t) with (e(t)e(tO) =2D5(t-t'). 



(21) 



Firstly let us solve the deterministic equation of motion under the influence of the different 
potentials depicted in Fig. HI All potentials offer their (global) minimum at x = 0, so this 
is the fixed point of the equation 

dx{t) 



dt 



For the Delta-starting distribution one can calculate the deterministic trajectory exactly 
x{t; jj,D,xo) = Q (t ^ 



VJTD 



Xq — sign(xo) \/ ^Dt 



There is a linear decrease (increase) of x to the fixed point x = in the time interval 
< t < l^ol/v^/lD. After reaching the fixed point the particle stays there if no noise is present. 

For the Gaussian distribution the drift term is highly nonlinear and the deterministic 
equation is not solvable exactly. In Fig. El the direction field of the solution of the deter- 
ministic motion is illustrated for different starting values, where the numerical calculation is 
done with the Runge-Kutta 4*^^ order algorithm with a step length h = 0.05. As seen from 
Fig. 131(b) the potential exhibits two different regimes, namely a parabolic one if |x| <C 1 and 
linear regime for |x| ^ 1. Assuming a piecewise approximation for the drift term as shown 
in Fig. ini the trajectories could be calculated exactly. Both pieces are merged together at 
the new parameter xq 



Xa{t) 




for t < ti 
for t > ti 

for t <t2 
for t> ti 



if Xo < 



//TP 



if 



/JTD 



< Xo < 

if xo>-4^ 



fJTD 




//////////////////// 
//////////////////// 
//////////////////// 
//A///////////////// 
//////////////////// 
1 n ! n 11 1 n I n 11 11 ! I 
1111111111111111111 



FIG. 5: Numerical calculated trajectory for "p^^ix) = exp(— x^), D 
-2;-l;0;l;2 



1 and xo 




FIG. 6: Comparison between the drift term /(x) (solid) and the piecewise approximative drift 
term /a(x) (dashdotted) forA = l,/i = l,Z) = l 

where t\ = — ^■^jy/pTD and ^2 = + ^'^/^/JTD. In the outer area the decrease of x{t) is 
linear, whereas after the cross-over it is exponential in the inner area. Because the piecewise 
approximation is quite rough we have studied a better one by replacing the drift term by 



f{x) = \/ [iD tanh 



X 



(22) 



where both regimes are matched, the long distance regime x ±oo as well as the short 
distance regime x — 0. The approximative drift term leads to the following solution for 
the trajectories 



arsinh 



sinh 



Xo 



which are plotted in comparison to the numerically calculated trajectories in Fig.|Hl Approx- 
imation (|^ underestimates the exact drift term that leads to an overestimated absolute 
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(a) f{x) and (b) Relative error 

faix) m 

FIG. 7: The drift term f{x) and the approximative one p2jl and the relative error (percentage) 
for A = 1, /X = 1 and D = 1 




FIG. 8: Trajectories x(t) (sohd) and the approximative trajectories Xa{t) (dashdotted) A = = 
= 1 and -4 < xo < 4 

value of x{t) in the intermediate regime, whereas for the piecewise approximation the drift 
term is overestimated and (the absolute value of) the deviation between x{t) and Xa{t) in 
the cross-over area is underestimated. A similar situation is observed for the exponential 
decreasing initial condition. Therefore, we skip this part. 

Concluding our paper we have discussed a simple evolution equation with a long-range mem- 
ory. Due to the permanent coupling to the initial distribution the system offers a stationary 
solution which depends on the special choice of the initial distribution. Then we have demon- 
strated that this evolution equation is fully related to a FPE with non-trivial deterministic 
forces. As a new feature it results that the deterministic part is also characterized by the 
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diffusion constant D, with other word the memory induces a stochastic behavior also within 
the deterministic part. As the consequence already the deterministic part of the underlying 
Langevin equation exhibits special trajectories. 
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